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Abstract 

By using a non-perturbative expansion and the dispersion relation for the Adler 
D-function we propose a new method for constructing the QCD effective coupling 
constant in the timelike region. 



The perturbative expansion is a powerful tool for performing calculations in quantum 
chromodynamics. Perturbation theory with renormalization group improvement has been 
widely applied to the description of various processes in both the spacelike and timelike 
regions. However, as is well known, the running coupling constant as(Q 2 ), the expansion 
parameter in QCD perturbation theory, is usually defined in the Euclidean range of mo- 
mentum. If one wishes to know the "running coupling constant" in the timelike region, 
and how it evolves, one must make use of some analytic continuation from the Euclidean 
( spacelike ) region to the timelike region. To this end one usually applies the dispersion 
relation for the Adler D-function [1] 

d s 



D{q 2 ) = ~q 2 r-*^R( 8 ), (1) 
Jo (s - a 2 ) 



where 



D(g 2 ) = <? 2 (- ^) n(<? 2 ) (2) 



and the function H(q 2 ) is defined by 

i J d 4 xexp(iq.x)(0\T{J^x)J u (0)}\0) oc (q^q u - g^q 2 )U(q 2 ) (3) 



and 

R{s) = — [n(s + ie) - n(s-ie)] = -Imn(s) 
2tt[ tx 
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The assumptions for writing these equations are the following: the function Il(g 2 ) is an 
analytic function in the complex g 2 -plane with a cut along the positive real axis, and the 
dispersion relation for the D-function does not require any subtractions. The contour in 
Eq. (4) goes from the point z = s — i e to the point z = s + i e and lies in the region 
where the function D(z) is an analytic function of z ( see Figure 1 ). 
Let us represent the functions D(q 2 ) and R(q 2 ) in the form 

D{q 2 ) oc J2Q 2 f n + do\t(q 2 )+d 1 \ 2 (q 2 )+...] 
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and 

R(s) oc £Q?[1 + n>A s ( S ) + nX 2 s (s) +■••] 
/ 

ee E^/ti + ^oAf^)], (6) 
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where we use the definition A = 0:5/(477). The index t means "t-channel" ( the spacelike 
region ) and s means "s-channel" ( the timelike region ). The coefficients d in Eq. (5) 
and r in Eq. (6) are numerically equal; thus Eq. (1) leads to a dispersion relation relating 
the effective coupling constants in the t and s channels: 

A?V) = -S 2 /°°-^Af( S ). (7) 
J (s - q 2 ) 

The solution of Eq. (7), given the analytic properties of A£ ff (t), is 

Af( S ) = -^ / S+ie -Af(.) (8) 

Z7T1 Js-ie Z 

Thus, Eq. (8) serves to define the running coupling constant in the timelike region. If 
we use the leading one-loop approximation for the running coupling constant A£ ff in the 
i-channel, 

where 60 = 11 — 2n//3 is the first coefficient of the /3-function and rif is the number of 
flavours, we obtain 

1 71 1 1 7T 2 

&o Af (s) = — arctan — ; — — — ■ = - — ; — — — - — g- — — — +••• (10) 

s V ; 7T ln(s/A 2 ) ln(s/A 2 ) 3 ln 3 (s/A 2 ) V ; 

Thus, one might conclude that the main effect of the analytic continuation from the 
spacelike to the timelike region is an additional contribution ( the so-called 7r 2 -term ) 
of order l/ln 3 (s/A 2 ). The 7r 2 -terms play an important role in the analysis of various 
processes. For example, in the r-lepton region, where 0:5 ~ 0.35 ( see, for exam- 
ple, a discussion in Ref. [2] ), the 7r 2 -correction is about twenty percent and the ratio 
R(s) = X s (s) / X t (q 2 ) ~ 0.8 for s = - q 2 = M 2 . The 7r 2 -terms and the problem of the 
construction of the running coupling constant in the timelike domain have been discussed 
in detail in Refs. [3-7]. 

However, the one-loop approximation to Af ff breaks the analytical properties of Af ff (z) 
for which we wrote the solution (8). Indeed, the function \^~ loop has an infrared singularity 



at z = — A 2 , which contradicts the assumption that the function \t(z) has only a cut in 
the complex z-plane. A consequence of this is the fact that Eq. (7) does not reproduce the 
original one-loop coupling constant in the t-channel if we use for A^ ff (s) the "solution" in 
the form of Eq. (8). In order to maintain the analytic properties of A cff one must make 
some modification to the initial coupling constant, which is to be substituted into Eq. (8). 

From a phenomenological point of view such modifications can be made in an ad hoc 
manner, with no fundamental justification. In the present note, however, we construct 
the running coupling constant in the timelike domain in the well- motivated framework 
of the non-perturbative method proposed in Ref. [8], which allows one systematically to 
determine the low energy structure in QCD. The method is based on the introduction 
of a new small expansion parameter. For small coupling constant cts{Q 2 ) the series thus 
obtained coincides with standard perturbation theory, so that high energy physics is 
essentially unaffected. Moreover the method continues to be valid in the nonperturbative 
region, where as(Q 2 ) > 1, since the expansion parameter in this approach remains small. 

The approach is related to the method of the Gaussian effective potential [9-12], the 
linear ^-expansion [13,14] and the approach based on variational perturbation theory 
[15,16]. Within these approaches the investigated quantity is written from the very be- 
ginning in the form of a series which provides a well-determined algorithm for calculating 
corrections up to any order. In addition, the series contains some auxiliary parameters 
which allow one to control the convergence properties of the series and to construct an 
expansion based on a new small parameter, as in this paper. The exact value of a given 
quantity is independent of these parameters. However, the corresponding approximate 
result contains a residual dependence upon those variables, and we may choose the latter 
so as to provide the best approximation to the quantity. A similar situation holds in stan- 
dard perturbation theory for the total cross section for e + e~ annihilation into hadrons: 
the physical quantity R e + e - is independent of a renormalization scheme, whereas a fi- 
nite perturbative approximation to R e + e - is scheme dependent ( in this connection see 
Refs. [17,18] ). 

The method described in Ref. [8] starts with the standard action of QCD, written as 
S(A,q,<p) = S 2 (A) + S 2 (q) + S 2 (y?) + gS 3 (A,q,(p) + g 2 S 4 (A) , (11) 
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where 82(A), including the gauge-fixing term, S 2 (q) and S 2 (ip) are the standard free 
actions of the gluon, quark and ghost fields, while S 3 (A, q, ip) contains the trilinear vertices 
and S±(A) the four-point gluon vertex. By use of an auxiliary field x»v the latter action 
can be rewritten as a trilinear action between A and x, so that S becomes 

S(A,q,p, x ) = S 2 (q) + S 2 (ip) + S 2 ( X ) + S(A, X ) + gS 3 (A,q,p). (12) 

We then introduce a new split between the free and interaction parts of the action 
according to 

S(A, q, p, X ) = S' (A, q, p, X ) + S' T (A, q, ip, X ) , (13) 

where 

S' = C 1 [S(A,x)+S 2 (q)+S 2 (p)]+r 1 S 2 (x), (14) 

and 

S'j = gS 3 (A,q, V )-(C 1 -l)[S(A, X ) + S 2 (q) + S 2 (p)} - (C 1 ~ 1) S 2 ( X ) , (15) 

thus introducing two variational parameters £, (. However, it turns out [8] that the two 
parameters must be related by £ = ( 3 in order to preserve gauge invariance. After a 
rescaling of the fields we obtain the following expression for a general Green function: 

G(- ■■) = j D QCD (■ ■ ■) V(A, q, ip) exp(* S ) , (16) 

with 

n I o n-k -k -1 

v(A, q ,ip) = EEj-^i-^) -n 



^ (n-k)\ K 8k' k\ [1 + k (C _1 - 1) ] u/2 
x [ 93 S 3 (A,q,ip)] k exp{t[glS 4 (A)]}, (17) 



where 



93 = [i + Mc 9 '-Dr • 94 = + -i)] i/2 (18) 

Here S'0 is the standard free action of QCD and k is a parameter introduced for convenience 
which is set equal to 1 at the end of the calculation. 

Analysis of the structure of this variational perturbation series shows that it can be 
organized in powers of the new small parameter a = 1 — £ if the standard coupling constant 
g is related to a by 

g 2 1 a 2 



where C is a positive constant. It is clear that for all values of A > the expansion 
parameter a obeys the inequality < a < 1. The variational parameter C is fixed, for 
example, by requiring the correct properties of the /3-function in the infrared region, as 
is discussed in more detail below. 

The result of the expansion for the expression (17) to order a 3 is 



Here we must make some comments. The one-loop level that we use here allows us to 
calculate the renormalization constant Z\ both to 0(a 2 ) and to 0(a 3 ). The expression 
(21) gives the 0(a 3 ) result, while P^(X) is obtained by simply omitting the term 9 a in the 
final bracket of Eq. (21). We will apply the phenomenon of induced convergence whereby 
the variational parameter is not fixed for all orders but rather depends on the order of 
the expansion. As was empirically noticed in Ref. [19] the results seem to converge if one 
changes the variational parameter from order to order using some variational principle. 
The induced convergence phenomenon has been discussed in detail in Ref. [20]. In Ref. [14] 
the convergence of the optimized (^-expansion has been proved in the cases of zero and 
one dimensions. 

The constant C is determined by insisting that at an appropriately large value of A the 
/^-function should behave like — /5(A) /A ~ 1, which ensures the correct infrared behaviour 
of the running coupling constant: as(Q 2 ) ~ Q~ 2 ■ This singular infrared behaviour of 
the coupling constant leads to the linear part of the quark-antiquark potential which is 
in agreement with the phenomenology of hadron spectroscopy. In 0(a 2 ) and 0(a 3 ) we 
find = 0.977 and = 4.1 respectively [21]. 

By using the renormalization group equation 




where A 3 = (iS 3 )/P, A 4 = (iS 4 )/P 2 and P = C/(4n) 2 . 
Within this approach the one-loop /3-function is 




2 (2 + 9a). 



(21) 




(22) 
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we obtain 



91 

11 Ql 2b 



/W(a) - / (i) (a ) , (23) 



where 



\ 2 12 B1l l-a 9 
f 2 \a) = - + — + 21 In. 



a 2 a a 1 — a 

rfiw \ 2 6 „ n , 18 1 624, „ , 5184, „ 9 , 

/ (3) («) = ^--- 481na -TTr^ + ^T ln(1 " a) + i^r ln(1 + 9 a) - (24) 

a a 11 1 — a 121 121 2 

The parameter ao and the momentum Qo in Eq. (23) are defined by some normalization 
point for our effective coupling constant 

A?V) = ^(l + 3a). (25) 

The effective coupling constant \f\q 2 ) is obtained from Eq. (25) by omitting the term 3a 
in the final bracket and using the constant instead of C^ 3 \ To fix a we use the value 
of the running coupling constant measured at the r-lepton mass [2]: Q = M T = 1.777 
GeV and as(M T ) = 0.35. In Figure 2 we show the behaviour of our running expansion 
parameter a(Q 2 ) in the i-channel, obtained from (23) using f^(a). As one can see, this 
parameter is indeed small in a wide range of Q 2 , and approaches 1 for very small Q 2 . 
Moreover, the effective expansion parameter in the method is actually a 2 (Q 2 ), with the 
0(a 4 ) coming in at two-loop level. Thus, the real expansion parameter is more or less 
that shown in Figure 2. 

Note that at large Q 2 = — q 2 the effective coupling constant (25) essentially coincides 
with the perturbative coupling constant. However, at small Q 2 their behaviours are 
quite different. While the perturbative coupling constant has a singularity at Q 2 = A 2 , 
the effective coupling (25) is finite at all values of Q 2 > 0. As follows from Eq. (23), 
the effective coupling constant (25) has a cut along the positive g 2 -axis. Thus, our 
approximation for the effective coupling constant in the t-channel does not contradict 
the dispersion relation (7), and we can make use of Eq. (8) to find the effective coupling 
constant in the s-channel. 

By using Eq. (8) one finds 

a?(*) = ^iak^W - tf (, V)], (26) 



7 



where 

(2) (a) = j-^—iZ ~ Ha - 4(1 -a) In a + 3(1 -a) In (1- a)] (27) 

and 

\ „ , 72 1 318 , „ 256 , „ 9 , 

^> (a ) = -4 In „ - -— + — - „) + — ln(l + -a) . (28) 

The values of a± obey the following equation (from Eq. (23)) 

f(a ± ) = f(a ) + 2 -^[\n^±in]. (29) 

In order to check our calculations we substituted X s (s) obtained from (26) into Eq. (7), 
and found good agreement with the initial Af ff (g 2 ). In Figure 3 we show the behaviour 
of Rea(s) and — Imc(s) in the s-channel, working to 0(a 3 ). For 0(a 2 ) we obtained a 
similar result. In Figure 4 we show the behaviour of the ratio R(s) = \^ (Q 2 ) / \[ l \q 2 ) 
with s = Q 2 = — q 2 versus s for % — 2 (dashed line) and % — 3 (solid line). The 
value of R(s) for s = M 2 is about 0.93, which represents a smaller correction than 
the perturbative 7r 2 -correction of Eq. (10). The agreement between successive orders is 
confirmation of the mechanism of induced convergence. As follows from the results we 
have obtained, the effects of analytic continuation can be essential for the definition of 
the QCD coupling constant from timelike processes, for instance from r-decay. We plan 
to discuss this problem in future publications. 
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Figure Captions 
Figure 1 The integration contour in Eq. (4). 

Figure 2 The behaviour of our running expansion parameter a(Q 2 ) in the t-channel. 

Figure 3 The behaviour of Rea(s) and — Ima(s) in the s-channel at 0(a 3 ). 

Figure 4 The behaviour of the ratio R(s) = A^(Q 2 )/A^(g 2 ) with s = Q 2 = — q 2 
versus s for % = 2 ( dashed line ) and % — 3 ( solid line ). 
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